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A 10,000 YEARS OLD EXPLOSION IN DR21 
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ABSTRACT 

Sensitive high angular resolution (~ 2") CO(2-l) line observations made with the Submillimeter Array 
(SMA) of the flow emanating from the high-mass star forming region DR21 located in the Cygnus X molecular 
cloud are presented. These new interferometric observations indicate that this well known enigmatic outflow 
appears to have been produced by an explosive event that took place about 10,000 years ago, and that might be 
related with the disintegration of a massive stellar system, as the one that occurred in Orion BN/KL 500 years 
ago, but about 20 times more energetic. This result therefore argues in favor of the idea that the disintegration 
of young stellar systems perhaps is a frequent phenomenon present during the formation of the massive stars. 
However, many more theoretical and observational studies are still needed to confirm our hypothesis. 
Subject headings: stars: pre-main sequence - ISM; jets and outflows - individual: (DR21) - individual: (DR21 
Outflow) 



1. EVITRODUCTION 

Located in the Cygnus X complex at a distance of 1.36 + 
0.12 kpc (Rygletal. 2012), the outflow in DR21 is proba- 
bly one of the most enigmatic high-velocity molecular out- 
flows associated with star forming regions. This outflow is 
extremely massive (> 3000 Mo) and energetic (> 10'*'* erg) 
with a luminosi ty in the 2 p^m band alone calculated to be 
1800 Lo (GardeiietaITT99Tt lGarden & C arlstrom"1992'). The 
DR21 outflow has been suggested to be energized by a mas- 
sive protostar with bolometric luminosity of about 10^"^ Lq 
(iGarden et al.ll991l:lGarden & Carlstromll992.) . which would 
exc eed the total bolom etric luminosity of DR21 region (10^ 
Lrai lHarvev et al.|[l977h . Howev er, such a massive star has so 
far never been clearly found out (jCruz^jonzalez et al. 2007). 

The outflow appears to be bipolar with its strongest molecu- 
lar lobes extending northeast-southwest, and emanating from 
a dense dusty filame nt that extends from north to south 
dSchneider et al.l 12010 1). However, there are some faint H2 
filaments that show different orientations with some of them 
pointing back to the DR21 main. 

It has been suggested that the axis of this bipolar out- 
flow might be located very close to the plane of the sky 
(e.g. i Garden et al.lll991l) . however, the high velocity wings 
obsei-ved in H2 and CO (-80.82 to +46.84 km s ') and 
the ov erlapping of blue/redshifted emission on both lobes 
dCruz- Gonzalez et al. 2007; Schneider et al. 2010) indicate a 
different scenario for its orientation. Commonly bipolar out- 
flows with their axes located close to the plane of the sky 
show slow radial v elocities with the ir blue/redshifted lobes 
well separated (e.g. i Pech et aT]|2012 b. However, overlapping 
blue/redshifted emission on both lobes might result from an 
outflow in the plane of the sky (that seems to be unlikely due 
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to its high velocity wings) if the outflow is thought of as a cone 
with one side of the cone coming toward us and the other side 
moving aw ay. 

Garden et al.l 
Cvgan owski et al.l 



fT99l : IGarden & CarlstromI (IT992h : 
(I2003h : iDavisetalJ (1200 71) additionallv 
reported the presence of some collimated flows emanating 
from DR21 with diflrerent orientations to the main NE-SW 
DR21 outflow. 

In high-density protostar clusters, envelopes and disks pro- 
vide a viscous medium that can dissipate the kinetic energy 
of passing stars, greatly enhancing the probability of cap- 
ture. This process opens the possibility of having protostel- 
lar mergers that may generate impulsive wide-angle outflows, 
shock-induced masers, radio continuum emission, and run- 
away massive protostars. Both the ejection of massive stars 
and the launch of the impulsive outflow may be the result 
of the dynamic interaction and rearrangement of a system of 
massive stars. T his process has been extensiv ely discussed in 
iBallv & Zinneck er (2005): lBallv et all (120111) . 

Very recent sensitive '^CO(2-l) millimeter and infrared ob- 
servations have suggested that the massive (10 Mq) and ener- 
getic (~ lO'*^ erg) outflow located in the Orion BN/KL region 
appears to have been produced by such a violent explosion 
during the disruption of the massive young stellar system of 
which the infrared and radio s ources BN, n, and / were a ll 
members about 5 00 years ago (iZaoata et al.ll2009L IMTallbl: 
iBallv et £111201 1[) . This suggests that the complex of CO and 
H2 emission in BN/KL is due to a phenomenon different from 
the well known collimated outflows produced by the young 
stars commonly associated with star formation. With this re- 
cent result in mind, one could consider if the flow emerging 
from DR21 could be a similar case as that occurred in BN/KL 
500 years ago. 

We thus here report new SMA interferometric observations 
of this outflow trying to elucidate its nature. 



2. OBSERVATIONS 
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Fig. 1. — High velocity redshifted (red/brown) and blueshifted (blue) CO(2-l) filaments in the DR21 outflow as obsei"ved with the SMA, overlaid on an Ht image from Davis et 
al. (2007). The blue contours are from 13% to 90% in steps of 3%, the peak of the molecular emission; the peak emission is 21 Jy beam"' km s"' . The red contours are from 12% to 
90% in steps of 5%, the peak of the molecular emission; the peak emission is 16 Jy beam"' km s"' . Each line represents one sequence of positions at which the CO emission peaks in 
consistent velocity channels. Most of the filaments point approximately toward the same central position (represented here by a box) q-=20'' 39"' 1 . 1 ^ ± 5" and (5=+42° 19' 37.9" ± 5" . 
The size of the box represents the eiTors given in the text. There are also some H2 faint filaments that appear to point toward the same central position IDavis et aU2007l) . These H2 
filaments are marked with white circles. We used the mosaicing mode to cover the entire DR21 outflow as far as it has been mapped in H2. The synthesized beam is shown in the lower 
left corner. 

The observations were made with the Submillimeter Array ^ 
(SMA) during 201 1 August and 2012 July/August. The SMA 
was in its subcompact, compact, and extended configurations 
with baselines ranging in projected length from 7 to 160 m. 
We used the mosaicing mode with half -power point spacing 
between field ce nters and co v ered the entire DR21 outflow 
mapped in H2 by iDavis et aTl (12007 ). see Figure 1 . We con- 
catenated the three data sets using the task in MIRIAD called 
"uvcat". The three diff'erent observations were identical, and 
only the antenna configuration of the SMA changed. 

The SMA correlator was configured in 24 spectral "chunks" 
(or windows) of 104 MHz each, with 128 channels distributed 
over each spectral window, providing a spectral resolution of 
0.8125 MHz (1.05 km s"') per channel. The receivers were 
tuned to a frequency of 230.5387 GHz in the upper side- 
band (USB), while the lower sideband (LSB) was centered 
on 220.5387 GHz. The '^€0(2-1) transition was detected in 
the USB at frequencies close to a frequency of 230.5 GHz. 
The full bandwidth of the SMA correlator is 8 GHz (4 GHz in 
each band). 

The zenith opacity {TriocHz) was ~ 0. 1 - 0.3, indicating rea- 
sonable weather conditions. Observations of Uranus provided 
the absolute scale for the flux density calibration. Phase and 
amplitude calibrators were the star MWC349a, and the quasar 
J2007H-404. Further technical descriptio ns of the SMA and its 
calibration schemes can be found in Ho et al. (2004). 

The data were calibrated using the IDL superset MIR, orig- 
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Astrophysics, and is funded by the Smithsonian Institution and the Academia 
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inally developed for the Owens Valley Radio Observatory 
( IScovilleetal.lll993h and adapted for the SMA.^ The caH- 
brated data were imaged and analyzed in a standard manner 
using the MlRlAD, and KARMA (Gooch 1996) packages. 
The line image rms noise was around 100 mJy beam"' for 
each velocity channel at an angular resolution of 2!'17 x 1!'89 
with a RA. = +74.0°. 

3. RESULTS AND DISCUSSION 

In Figure 1, we show the '^CO(2-l) molecular emission 
detected in our SMA mosaic towards the DR21 region. We 
construct this large scale integrated intensity map only us- 
ing the emission arising from the high velocity gas and thus 
avoiding the emission in the range from -15 to 15 km s"'. 
In this velocity range the line emission is spatially extended, 
and could not be well sampled by the SMA. In our channel 
velocity maps, we found that the redshifted '^CO(2-l) emis- 
sion shows radial velocities from -1-3 up to +50 km s"', while 
the blueshifted one from +3 down to -80 km s This ve- 
locity ran ge is in very good agreement to that p reviously ob- 
served by iGarfenitdMllir; ISchne^^ How- 
ever, our high angular resolution observations revealed clear 
molecular filaments with different orientations, and some of 
them well correlated with H2 emission (Figure 1). 

The CO filaments reported here are composed of about 
three hundred compact features observed in our channel ve- 
locity cube with velocity windows of 1 km s"' and their po- 
sitions on the sky are shown in our Figure 1. The observa- 
tional parameters of these compact features were found us- 



^ The MIR cookbook by C. Qi 
|http://cfa-www.harvard.edu/~cqi/inircook.htmI 
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Fig. 2. — Position-velocity diagram of the jet-like CO complex: Radial velocity as function of on-the-sky distance from the common center for each of the 25 
C0(2-l) filaments, with blueshifted structures shown in blue, and the redshifted ones in red. Nearly all filaments seem to start from the same common radial 
velocity of about -3 km s"' (the small black circle). Black gradient dashed-lines start from -3.0 km s"' . Note that velocities between -15 and +15 km s"' 
could not be studied because the interferometer cannot sample well the extended molecular gas. 



ing the task imfit of MIRIAD. Each of these filaments show 
clear velocity gradients that follow a Hubble law (Figure 2). 
Furthermore, most of these filaments appear to point out a 
common velocity center located at about -3.0 km s"', in very 
good agreement with the ambient molecular velocity of DR2 1 
(-3.0 km s-': ISchneider et al.ll2010l) . 

Many of these physical features are also found in the molec- 
ular outflow located in BN/KL and are related to an explosive 
event occurred 500 years ago as described in the introduction. 
However, the outflow in DR21 has a clear east-west bipolar 
shape observed on the H2 maps. This H2 bipolar shape could 
be due to the highly extinction produced by the north-south 
dusty lane where many of the CO filaments are revealed for 
the first time. This al so seems to be the c ase in the outflow 
located in BN/KL, see lZapata etalj ( 1201 lal) . 

Using a similar approach to that applied in 'Zapat a et al.l 
(^009), we determine the position of the origin of the CO fil- 
aments in Q'(J2000)=20'' 39"' 1.1' + 0.3' and 5(2000)=H-42° 
19' 37.9" + 5". 

This position coincides well with the center of an ex- 
panding cometary HII region (see Figure 3) imaged by 
ICyganowski et al.l (120031) . This positional coincidence sug- 
gests that both the explosive outflow and the large cometary 
HII region were probably produced by the same mechanism. 
Moreover, at the position of the outflow origin there is not a 
single young (radio, submillimeter, or infrared) star, suggest- 
ing that maybe its "source" is no longer located the re, as in the 
case of the explosive outflow in Orion BN/KL ( Zapata et in 
|2009). There is also a depression of centimeter continuum 
emission at this position (see Figure 3). 

The explosive Becklin-Neugebauer (BN)/Kleinman-Low 
(KL) outflow emerging from OMCl behind the Orion Neb- 
ula was powered by the dynamical decay of a non-hierarchical 
multiple system ~ 500 years ago that ejected the massive stars 
I, BN, and source n, with velocities of about 10-30 km s ' 



(Zapa ta etal.ll2009l 12011 alibi |Bally et al."2011'). In this dy- 
namical decay of a non-hierarchical multiple system a kinetic 
energy of about 10^' ergs was liberated. 

If we assume that the outflow in DR21 was produced by 
a dynamical decay of a non-hierarchical multiple system, the 
kinetic energy liberated on that disruption was about 20 times 
larger than that in t he Orion BN/KL system (> 2 x 10"*** erg; 
IGarden et al.l[T99Tt IGarden & Carlstromll 19921) . 

Again, very crudely assuming that the velocities of the stars 
ejected in DR21 have similar velocities as the ones located 
in Orion BN/KL (this without any evidence), and a dynami- 
cal age of about 10,000 years'* (obtained from our molecular 
data, tj ~ 1.2 X 10'^ cm / 4 x 10'' cm s"' ~ 10, 000 yrs) for 
the outflow in DR21. The reader should note that this estima- 
tion is also an approximation, since the values for the velocity 
are only radial velocities, and the displacements are only the 
projected components. 

We thus found that the ejected stars should be located about 
25" far from the center of the DR21 outflow. This is exactly 
the distance of the bright infrared source "DR21-D" associ- 
ated with a small cometary HII region (Figure 3) with re- 
spect to the origin of the outflow in DR21. This suggests that 
maybe DR21-D is related with the explosive event in DR21. 
The Figure 3 also reveals a reflecting east-west cavity associ- 
ated with the large cometary HII region. These infrared neb- 
ulosities are also observed in Orion BN/KL and most of them 
are not heated by stars. They are possibly heated by strong 
shocks produced by the explosive outflow dMenten & ReidI 
IB95; Zapata et al. 2011a). 

All the characteristics found in these observations in the 
outflow located in DR21 are very similar to the ones re- 
ported in the expl osive outflow located in Orion BN/KL 
(iZapata et al.ll2009l) . We thus propose that the complex of 

This is an approximate 1 st order estimation as in the case of the Orion 
BN/KL iZapata et al..2009.) 
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Fig. 3. — Redshifted and blueshifted CO filaments in the DR21 outflow as observed with the SMA. overlaid on a Spitzer infrared color image and a VLA 6 cm. emission contour 
image. In the Spitzer infrared color image composite the red represents the 8.0 fim, green the 5.8 fim, and blue the 3.6 fim emission. The contours are 5% to 90% with steps of 10% the 
peak of the centimeter emission; The peak of the centimeter image is 0.13 Jy Beam"^ . The center of the explosive molecular outflow is marked with a white square. Note that this center 
coincides with the center of the compact HII region revealed by the 6.0 cm. continuum emission map. Additionally, the infrared map reveals a continuum source associated with the 
cometary HU region called DR21-D. Note that the white square does not lie at the peak of the 6.0 cm. emission, rather it is located in the middle of the HII region, where no emission 
is detected at all. 



CO and H2 emission located in DR21 might have been pro- 
duced by an explosive phenomenon different from the stan- 
dard disk-outflows associated with star formation. This out- 
flow was probably originated during the disintegration of a 
massive stellar system as the one occurred in Orion BN-KL 
some 500 years ago. As mentioned in the Introduction sec- 
tion, the process of that disintegr ation of the stellar system 
has been discussed in detail by Bally & ZinneckeJ (120051) . 
The presence of explosive phenomena in two of the best stud- 
ied regions of massive star formation suggests that the phe- 
nomenon may be more common than previously thought and 
that future detailed studies may reveal more of them in other 
regions. 

4. CONCLUSIONS 

Sensitive high angular resolution (~ 2") CO(2-l) line ob- 
servations made with the Submillimeter Array of the flow 
emanating from the massive star forming region DR21 lo- 
cated in the Cygnus X molecular cloud have been presented. 
We found about 25 molecular blueshifted and redshifted fila- 



ments. These molecular filaments follow nearly straight lines 
and appear to all point toward a common center The radial 
velocity along each filament follows clearly a Hubble law, 
indicative of an explosive event. Our observations suggest 
that this outflow appears to have been produced by a violent 
explosion that took place about 10,000 years ago, and that 
seems to be related with the disintegration of a massive non- 
hierarchical stellar system, as the one that occurred in Orion 
BN/KL, but about 20 times more energetic. However, many 
more theoretical and observational studies are still needed to 
confirm our hypothesis. 
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